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. Abstract 

O 

, One way of avoiding the complex action problem in lattice QCD at non-zero 

q \ density is to simulate QCD-like theories with a real action, such as two-colour 

QCD. The symmetries of two-colour QCD with quarks in the fundamental and 
Q ■ in the adjoint representation are described, and the status of lattice simula- 

tions is reviewed, with particular emphasis on comparison with predictions 
from chiral perturbation theory. Finally, we discuss how the lessons from 

J> \ two-colour QCD may be carried over to physical QCD. 
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I. INTRODUCTION 



Recently, there has been a considerable interest in QCD at non-zero chemical potential, 
after a number of model studies have indicated a rich phase structure [T|-|3| (for a review, 
see [§]). Clearly, it would be desirable if these predictions could be tested by first-principles, 
non-perturbative studies, e.g. lattice QCD. Unfortunately, lattice simulations of QCD at non- 
zero baryon density using standard methods are in practice impossible because the action 
(and the fermion determinant) becomes complex once the chemical potential is introduced, 
causing importance sampling to fail. It is possible to split the determinant into a modulus 
and a phase, simulating with the modulus of the determinant as the measure and reweighting 
■ the observables with the phase, 

_ ((Oarg(detM))) 
{ ' ((arg(detM))) ' 1 ' 

where ((■■■)) denotes the expectation value with respect to the positive real measure. How- 
ever, the denominator in ([T]) is effectively the ratio of the partition functions of two differ- 
ent theories: the true theory and one with a positive real measure. This should scale as 
exp(— AF), where AF is the difference in free energy between the two theories. Since the 
free energy is an extensive quantity, the computational effort required to obtain a reliable 
sample rises exponentially with the volume. 

A number of approaches have been tried to overcome this problem. In the Hamiltonian 
formalism, the problem does not arise. Analytical results have been obtained in the strong 
coupling limit [|5]||, but so far no method for numerical simulations exists. 

With an imaginary chemical potential JF], the action becomes real and positive, so sim- 
ulations are straightforward. The problem is whether an analytical continuation to real fi is 
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possible. It works at high temperature |8|||, where also other approaches may be success- 
fully employed 0, but does not seem to be possible at zero or low temperatures. Imaginary 
chemical potential can also be used to formulate a quenched limit of QCD in the background 
of a non-zero number of static quarks [T(j . 

Cluster algorithms JTT] may provide a way of eliminating the sign problem by summing 
analytically over configurations in a cluster in such a way that the contribution to the 
partition function from each cluster is always positive definite. So far, these methods have 
been applied to a number of spin models; however, an application to QCD has yet to be 
found. 

Finally, the sign problem may be avoided by simulating theories which resemble QCD, 
but have a real action even at non-zero chemical potential. One such theory is QCD at 



non-zero isospin density p2]JT3[] , which is of intrinsic interest because it corresponds to part 
of the phase diagram for asymmetric nuclear matter. Another class of theories encompasses 
two-colour QCD with fermions in the fundamental representation, as well as QCD with 
adjoint fermions, for any number of colours. The remainder of this review will focus on 
what can be learnt from lattice simulations of these theories. 



II. THEORIES WITH REAL ACTION 



The chemical potential /i is introduced on the lattice by multiplying the forward timelike 
links by and the backward timelike links by e~ iX ]14]]. It can be shown jn|[nj that the 
determinant det M of the fermion matrix M in two-colour QCD is real, even for non-zero 
/i, both in the continuum and on the lattice. However, it is only possible to demonstrate 
that it is positive 
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in the cases of continuum or Wilson adjoint fermions and staggered 
fundamental fermions. Indeed, we will see in section [TV] that in the case of staggered adjoint 
fermions there are configurations with a negative determinant, leading to a sign problem at 
large /i. 



A. Symmetry breaking pattern 

In the chiral limit, the action for two-colour QCD with N flavours has a U(N)l <8> U(N)r 
symmetry, which for staggered fermions is manifest as independent U(N) symmetries for the 
even and odd sites. At /x = this enlarges to a U(2N) symmetry. This can be seen most 
easily by introducing new fields, 

Xe = (Xe, -X T eT 2 ) X Q = ( ) (2) 

for (staggered) fundamental quarks, and 

X e =(Xe,X T e) X a = ( g ) (3) 

for adjoint quarks. The action can then be written as 
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S=V^)fcx)( ' n )l J ,A.r)X (.r v) (1) 

X e {x) 

where x e denotes the even sites. In the continuum, the equivalent fields are 

fundamental: * = ^ ^ ^„ adjoint: * = ^ (5) 

which gives the continuum lagrangian 

C = itfa v (D v -nB v )V B u = (I °A5 I/0 - (6) 



The explicity chiral symmetry breaking term in the Wilson fermion action means that there is 
no equivalent enlarged symmetry for Wilson fermions; however, new fields may be introduced 
analogously to the continuum case, and the enlarged symmetries will be broken by 0(a) 
terms in the action. 

The chiral condensate can be written in terms of the new fields, 

(7) 

for staggered fermions, and 

W = * T °2%[ L "nM* + h -C- (8) 




in the continuum and for Wilson fermions. In both cases, the + sign is for fundamental 
fermions and the — sign for adjoint, while T is r 2 for fundamental fermions and 1 for 
adjoint. A nonzero chiral condensate thereby breaks down the U(2N) symmetry to 0(2N) 
for fundamental fermions and Sp(2iV) for adjoint fermions, giving rise to N(2N + 1) and 
N(2N — 1) Goldstone modes respectively. Of these, there will be iV 2 mesonic states, while 
the remaining iV(iV ± 1) will be diquarks. In the continuum, and for Wilson fermions, the 
pattern will be the opposite (modulo the 1 mode destroyed by the axial anomaly in the 
continuum), but for 1 flavour of fundamental quarks, there is no chiral symmetry in the first 
place. From this we see that in the case of iV = 1 adjoint staggered fermions, and only in 
this case, are there no diquark Goldstone modes. 

For m/0, all the pseudo- Goldstone modes remain degenerate, with masses oc \fm. 
As the chemical potential \i increases, the ground state will begin to be populated with 
baryonic matter. The transition to a ground state containing matter occurs when /i — /i ~ 
mb/n q , where nib is the mass of the lightest baryon, and this baryon contains n q quarks. At 
this point, the baryon number density n becomes non-zero, where n is given by 

n = I^e^ - l)U (x)i(>(x + 6) (9) 
+^j(x + 6)e^( 7o + l)Ui(x)ip(x)) 
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for Wilson fermions, and 



n = ^(x(x)vo(x)[e^Uo(x)x(x + 0) + e^U^x - 6) X (x - 6)]). (10) 

for staggered fermions. Where there are diquark Goldstone modes, those states will be the 
lightest baryons in the spectrum. This means that for most variants of two-colour QCD we 
expect /x D ~ r/ij/2, in contrast to the much larger value 171^/3 expected in real (three-colour) 
QCD. The exception is two-colour QCD with one flavour of adjoint staggered quarks. 

In the limit of small m and /x, the behaviour of (ipip) , the diquark condensate (ipip) , and 
n as functions of m and /x can be calculated in chiral perturbation theory. If we define the 
rescaled variables 

2/x ($ip) (-0-0) _ m n0 n 

x = ' y= 7T7V> Z = 7TJV^ n = g TTT\ ' ( n ) 

where the subscript denotes values at /x = 0, the prediction from x?T for the models with 
diquark Goldstone modes is 





* = {*h i\ ( 12 ) 



B. Diquark condensation 

At large chemical potential, the relevant degrees of freedom will be quarks with momenta 
near the Fermi surface. The attractive quark-quark interaction will give rise to instability 
with respect to condensation of diquark pairs at opposite sides of the Fermi surface. In 
physical QCD, the diquark condensate cannot be a colour singlet, so the gauge symmetry is 
spontaneously broken, giving rise to the phenomenon of colour superconductivity. 

In two-colour QCD, on the other hand, there may be the possibility of gauge singlet 
diquarks, which will be energetically favoured compared to non-singlet states. Indeed, in 
the previous section we saw that most variants of two-colour QCD have diquark Goldstone 
modes, which will be the preferred channel for diquark condensation. In the N = 1 staggered 
adjoint model, this is not the case, and we do not know a priori in which channel the 
condensation will occur. We must proceed by constructing possible operators which obey 
the Pauli principle and making additional assumptions about locality, Lorentz structure and 
gauge mvanance p| . At least one of the possible condensates constructed this way gives 



rise to a colour superconding ground state. 

The standard way of computing the diquark condensate on the lattice is to introduce a 
diquark source term into the action |T7|]. For two-colour QCD with fundamental staggered 
quarks the action then becomes 

Sf = J2 x(x)M xy x(y) [x T (x)r 2 x(x) + x(x)t 2 x t (x)\ (13) 

x,y x 

= (x,X T )(i?J*?)(?)=X T AU]X. (14) 
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M JT 2 
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In this case, the partition function becomes proportional to the Pfaffian PD4[j]. The diquark 
condensate (x Tr 2X) ma y De evaluated by taking 



An alternative approach |T8| is to rewrite the Pfaffian as 



PL4[j] = Pf(B+j) = ± v /det(B+j) (16) 

where 

This can be expanded as a polynomial in j by diagonalising B 2 , obviating the need to 
simulate at non-zero diquark source. Since this gives the Pfaffian at any j, it can also be 
used to determine the diquark condensate using the probability distribution function fT9|. 



III. SIMULATIONS WITH FUNDAMENTAL QUARKS 



In the past year and a half, a number of groups have been performing lattice simulations 
of two-colour QCD with fundamental staggered fermions both at zero |20| - |23| , |T3|1 and non- 
zero |2^j25|1 temperature. Also, one group is performing simulations with Wilson fermions 



Aloisio et al. pi, 22 have performed simulations in the strong coupling limit for a number 
of quark masses, flavours and lattice volumes. Fig. |I| shows results for the chiral condensate, 
the diquark condensate and the baryon number density, at m = 0.2 and a non-zero source 
j = 0.02, for two different lattice sizes. Also shown are the predictions from flT2]). The 
agreement between the prediction and the numerical results is quite striking, considering 
that this is far from the continuum limit. This suggests a weak f3 dependence. Fig. ^| shows 
the diquark condensate at zero diquark source, for Nf = 1 and a range of quark masses. 
Again, we see a very good agreement with the prediction ([12]) . At larger /i, we see that 
the value of {ipip) drops, going to zero at high \i. This second transition is due to lattice 
artefacts connected with the saturation of lattice sites with fermions. In the infinite volume 
limit it is expected to disappear. 

Simulations at non-zero /i with a diquark source have been performed by Kogut and 
Sinclair [O]. Results for the diquark condensate and chiral condensate are shown in fig. 
Again, we see the chiral condensate dropping and the diquark condensate rising for fi > 
in agreement with ;\PT. We also see the same large-/x saturation behaviour for the 
diquark condensate as in . Fig. [| shows results for pion and scalar diquark masses. The 
scalar diquark mass falls roughly as m n — 2fi as fi approaches m vr /2. The pion mass remains 
constant up to fi ~ rn n /2, after which it falls to zero. This is again in accordance with the 
expectation from %PT. 

The spectrum of the Dirac operator has been studied in some detail by the Vienna 
group @ for staggered fermions and by the Hiroshima group [2£| for Wilson fermions. A 
preliminary study of topology at non-zero temperature has also been performed 
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IV. SIMULATIONS WITH ADJOINT QUARKS 



As indicated in section |11 A| , two-colour QCD with one flavour of adjoint staggered 
fermions has features which makes it in some senses more 'QCD-like' than other variants of 
two-colour QCD. In particular, it has no diquark Goldstone modes, so we expect an onset 
transition at a value of the chemical potential different from m n /2 — possibly at /x = mjv/3 
where denotes the mass of the lightest three-quark baryon (the 'nucleon'). It also has a 
sign problem. 

The theory has been simulated [[U],^7]] usrn g two different algorithms: Hybrid Monte 



Carlo, which is not able to change the sign of the determinant, and therefore only simulates 



the positive determinant sector of the theory, and the two-step multibosonic algorithm p8 
which is able to take the sign properly into account. Figs |5| and |5] show the results from HMC 
simulations for y and n of ([□]) respectively, for a range of values for the quark mass m and 
chemical potential /x. Up to x ~ 1.5, the data collapse onto a universal curve, which agrees 
well with the predictions of flT2|). Even at larger x, the data for the chiral condensate lie 
close to the xPT prediction. At x 2, however, the data for different m diverge, indicating 



that %PT may be breaking down. Results for the plaquette p7| show a drop in its value for 
/i > m n /2, presumably due to Pauli blocking, while the pion mass appears to agree with 
the x?T prediction = 2/x at x > 2. 

The agreement between the predictions of %PT and these results paradoxically enough 
presents a problem, since this model is not supposed to contain any diquark Goldstone 
modes, and thus the %PT predictions of (|12"D are not valid in this case. In particular, 
there should not be any onset transition at /x — The suspicion must be that this 

contradiction is due to the fact that HMC does not change the sign of the determinant, and 
that it therefore simulates the wrong theory — a theory with conjugate quarks. 

The simulation points for the TSMB algorithm were selected to focus on the effect of the 
sign, with one point at /x = 0, one just past the HMC onset transition, and one deeper into 
the dense region. With this algorithm, a reweighting factor r and the sign of det M must 



be determined for each configuration [29|. The expectation value of an observable O is then 
determined by the ratio 

(O) = (QXrX ^ ) . (18) 
(r x sign) 

The results for (ijjijj) and n, together with the corresponding HMC results, are summarised 
in Table [J For TSMB at /i ^ we also include observables determined separately in each 
sign sector, defined by (0)± = (O x r)±J (r)±. At /x = 0.0 the two algorithms agree, as they 
should. Also, the results in the positive determinant sector for TSMB at larger /x agree with 
the HMC results. However, the results for the negative determinant sector are significantly 
different. This difference has the effect of bringing the average both for (ipip) an d for n back 
to values consistent with the /x = values. This is an indication that at /x = 0.36 and quite 
possibly also at /x = 0.4, the system is still in the vacuum phase, which means that the onset 
transition in this model occurs £it £l lcir ger \x than for other variants of two-colour QCD. 
This is consistent with the symmetry-based arguments of section [II A] that this model has 
no baryonic Goldstone modes. 
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V. CONCLUSIONS 



Substantial progress has been made recently in lattice simulations of two-colour QCD 
at non-zero density, with both fundamental and adjoint quarks. The simulations with fun- 
damental quarks nicely reproduce the predictions from chiral perturbation theory for the 
chiral condensate, diquark condensate, and baryon density, except at very large chemical 
potentials. Here, saturation effects are observed, especially for the diquark condensate. The 
meson and diquark spectrum is also being analysed, with preliminary results for the pion 
and scalar diquark masses again in rough agreement with chiral perturbation theory. 

Two-colour QCD with one flavour of adjoint staggered quark is not expected to have 
any diquark Goldstone modes, unlike all other variants of two-colour QCD. It also has a 
sign problem. Simulations of this model restricted to the sector with a positive fermion de- 
terminant reproduce the predictions of chiral perturbation theory for theories with diquark 
Goldstone modes, including an early onset transition at /i ~ m n /2 — although the break- 
down of xPT may be observed at larger /i. When configurations with negative determinants 
are included, we find a strong correlation between the sign and the value of observables. 
This effect appears to lead to a cancellation of the early onset transition. This observation 
may hold the key to understanding the problem of the premature onset which has bedevilled 
previous attempts to simulate physical QCD at non-zero density. 

The presence of diquark modes which are degenerate with the pion means that the 
study two-colour QCD is of limited usefulness when it comes to studying directly the onset 
transition, hadron spectrum and diquark condensation in physical QCD at non-zero density. 
However, useful experience may be gained by comparing the results of lattice simulations 
with those of other methods which are also applicable to physical QCD. Of particular interest 
would be the study of gluodynamics, where SU(2) and SU(3) are expected to exhibit similar 
behaviour, even at non-zero \i. Thus it might be possible to cast light on the deconfinement 
transition at high jj, and low T by studying two-colour QCD. 
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FIGURES 




FIG. 1. Baryon density, chiral condensate and diquark condensate vs. chemical potential, from 
pi[| . The solid lines are the predictions of (fL2|). 
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FIG. 2. Diquark condensate for a 6 4 lattice, Nf = 1, for m = 0.025 (diamonds), 0.05 (squares) 
and 0.2 (stars) at strong coupling, from [^2|, with the predictions of (|T^). 
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FIG. 3. Chiral condensate (left) and diquark condensate (right) vs. chemical potential, for one 
flavour fundamental staggered quarks, on an 8 lattice; from pi. The arrow indicates /i ~ 
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FIG. 4. Pion and scalar diquark masses as functions of chemical potential, for one flavour of 

TBtt — 2/i. 



fundamental staggered quarks, on an 8 lattice; from [13]. The straight line is m 
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FIG. 6. Baryon density vs. chemical potential for one flavour of adjoint staggered quarks, using 
the rescaled variables of eq. (|Tl|); from [27|. 
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TABLES 





n 




TSMB 




HMC 






(O) 


<o>+ 


<o>- 




(xx) 


0.0 


1.525(3) 






1.526(1) 




0.36 


1.551(10) 


1.521(8) 


1.176(37) 


1.485(9) 




0.4 


1.50(15) 


1.248(17) 


1.177(20) 


1.253(10) 


n 


0.0 


0.0000(13) 






-0.0002(3) 




0.36 


-0.0003(80) 


0.0199(64) 


0.252(32) 


0.0172(28) 




0.4 


0.07(10) 


0.177(14) 


0.208(15) 


0.1667(90) 



TABLE I. A comparison of results between TSMB and HMC for one flavour of adjoint stag- 
gered quarks. Due to long autocorrelation times, the errors in the HMC results are probably- 
underestimated. 
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